Reorientation dynamics in thin glassy films 
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We present a study of orientational relaxation dynamics in thin films of a low-molecular-weight 
glass-former as a function of temperature and film thickness. The relaxation is probed by second- 
harmonic generation after release of a poling electric field. From the measured decays of the second- 
harmonic signal and their fitting with a stretched exponential, we can determine the distribution 
of relaxation times in the system. As temperature decreases from above the glass transition, we 
observe that the width of the distribution first increases under confinement, but that deeper in the 
glassy state, confinement has no effect anymore on the dynamics. 

PACS numbers: 64.70.Pf, 42.65.An, 68.60.-p 



I. INTRODUCTION 

In the liquid phase of glass formers, the viscosity and 
any sort of relaxation time increase enormously upon 
cooling. Relaxation times eventually become longer than 
the experimental time scale, which is essentially the time 
during which the system is allowed to equilibrate after a 
change of an external parameter or constraint. As a con- 
sequence, no changes are observable anymore in the sys- 
tem and the structural configuration of the liquid appears 
to freeze below the so-called glass transition temperature 
T g . A deep understanding of the nature of the glass tran- 
sition it is still lacking and reaching it is a challenging 
goal. Many theories have been proposed 0, IS IS H 01 
and revisited, but all fail to explain the whole collection 
of typical features of the glass transition. Somehow each 
theory applies to a limited range of conditions, either 
thermodynamic or kinetic. 

One of the widely used theories of the glass transi- 
tion is that of Adam and Gibbs [6j, which is based on 
the concept of cooperatively rearranging regions (CRRs) . 
According to this model, as the temperature of a liquid 
decreases toward the glass transition, the motion of the 
molecules is no longer independent, but the dynamics 
becomes correlated over a certain number of molecules. 
The spatial extent of this correlation is generally defined 
as the cooperative length £. Parallel to the develop- 
ment of theories based on the idea of cooperative mo- 
tion, other models focused on explaining the fact that the 
dominant a-relaxation is highly non-exponential. This 
non-exponentiality can be related to the superposition of 
relaxation processes with different relaxation times asso- 
ciated with different regions in the material, the so-called 
dynamic heterogeneities. The concept is indeed related 
to the CRRs, which can be considered as fast relaxing re- 
gions surrounded by slow or immobile molecules 0. So 
the length scale of heterogeneities and the cooperative 
length should essentially be the same. Both simulation 
and experiments have been used to provide evidence for 
the existence of such a length scale and to evaluate it (see 
i e IS IS El El d) The obtained values at T g or just 
above it, are of the order of a few nanometers, both for 
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FIG. 1: Schematic representation of the effect of confinement 
on cooperative motion: as the size of the system decreases, 
cooperatively rearranging regions (the radius of which is in- 
dicated by the arrow) gets modified, hence also the dynamics 
of the system. 



low-molecular weight and polymeric systems. Theories 
predict that this length should increase as temperature 
decreases. 

The concept of cooperativity suggests that the dy- 
namics of glass-formers should be modified if the size 
of the system becomes comparable to the cooperative 
length £. In this size range, it is the size of the sys- 
tem that determines the size of cooperatively rearrang- 
ing regions (Fig. ^) and therefore the dynamic behavior 
of the system. So one expects modifications of the dy- 
namics, namely of the distribution of relaxation time and 
the glass transition temperature, in glassy systems with 
sizes in the mesoscopic range. Taking into account only 
pure geometrical effects (without interactions with walls) 
and assuming the presence of cooperatively rearranging 
regions of size £, we can expect the following changes 
in the dynamic behavior of glass-formers when they are 
confined in thin films ^3- For thicknesses rather large 
with respect to £, confinement only affects molecules in 
the vicinity of the surfaces (Fig. [I}. Surface molecules be- 
long to cooperatively rearranging regions that are smaller 
than the ones present in bulk. Therefore their relax- 
ation should be faster than that of molecules in the mid- 
dle of the system. As the size of the system decreases, 
the proportion of molecules affected by the surfaces in- 
creases. As a consequence, the distribution of relaxation 
times in the system should broaden toward the short time 
scales leading to a decrease of the mean relaxation time. 
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FIG. 2: Chemical structure of the glass- former. 



When the film thickness becomes comparable with £, all 
molecules are affected by confinement and the size of 
the cooperatively rearranging regions is determined by 
the size of the system. Decreasing the size even more 
down to a molecular size will make the system homo- 
geneous again, with all molecules having essentially the 
same environment. However, the relaxation time will be 
smaller than the value in bulk, because the motion of the 
molecules becomes essentially individual again. 

Research on the effect of confinement on glass dy- 
namics has concentrated on two types of systems: low- 
molecular-weight glass-formers in porous media and poly- 
mers in thin films. In low-molecular- weight glass-formers, 
as measured by dielectric spectroscopy, the characteris- 
tic time corresponding to the a orientational relaxation 
is generally decreasing and the distribution of relaxation 
times is broadening as the size of the system decreases 
(l~3. ITU llH Il7j . From the effects of size on the dynam- 
ics, it is possible to make an estimate of the cooperative 
length £. The value of £ for salol in the vicinity of the 
calorimetric glas s transition temperature was estimated 
to be > 70 A [15| . Similar values of £ at T g were obtained 
for propylene glycol, with £ < 60 A and for jV-methyl-e- 
caprolactam (£ < 120 A) |l8j|. Confinement in porous 
media, however, introduces curvature effects and is not 
well defined in size because of the polydispersity in pore 
sizes and the presence of pore junctions. These problems 
are avoided in thin films, which have so far only been used 
for polymeric systems. In free-standing films, the mean 
relaxation time was found to be significantly lower than 
in bulk while the relaxation dynamics followed the same 
stretched exponential as in bulk 19] . In contrast a sig- 
nificant broadening of the distribution of the relaxation 
times has been observed in films supported by a substrate 
pol l2l| . There has also been one study of the effect of 
confinement on the width of distribution of the relax- 
ation times as a function of temperature in a side-chain 
polymer poj ]. It was found that, at all temperatures, the 
average relaxation time is essentially independent of the 
film thickness, while the width of the distribution of the 
relaxation times increases as thickness decreases. The 
widening of the distribution was however found to be 
the same over the whole temperature range considered 
(T g - 13 °C < T < T g + 25 °C). Studies on polymers 
then have not shown so far any clear trend. Because 
of their chain character, polymers exhibits bulk dynamic 



processes and a behavior under confinement that is spe- 
cific to these materials rather than glasses in general. 

To overcome these problems, we have studied the dy- 
namic behavior of low-molecular-weight glass-formers in 
thin films. To our knowledge there is no study re- 
ported in literature on the dynamics of low-molecular- 
weight systems in as simple a confining geometry as thin 
films. We have chosen a hydrogen-bonded glass-former, 
with a fairly simple chemical structure (see Fig- El) j from 
which we can expect relatively simple relaxation pro- 
cesses. Moreover, small molecules can be obtained with 
a higher chemical purity and without molecular weight 
or composition distribution. The small size of our com- 
pound has however the disadvantage that the character- 
istic lengths in the system should be smaller than in poly- 
meric systems. As a consequence, confinement effects will 
appear at smaller film thicknesses. 

In previous measurements performed with x-ray reflec- 
tivity [2^, we have seen one effect of confinement on a 
glass former into a thin film: the glass transition signifi- 
cantly widens. This was interpreted as originating from 
an increase of the inhomogeneity of the dynamics in the 
system. To obtain a direct evidence for the variation of 
the inhomogeneity under confinement and study in de- 
tail how confinement affects the dynamics of the system, 
we have examined the reorientation dynamics of our low- 
molecular- weight glass- former in thin films, starting from 
an initial state imposed by an external field. The dynam- 
ics has been followed experimentally over a wide range of 
temperatures around T g , an extended time window and, 
for confining sizes down to a few molecular units, which 
is expected to be of the order of £. We have used opti- 
cal second-harmonic generation technique to follow the 
re-orientation process. Second-order non-linear optical 
effects, such as second-harmonic generation, are in prin- 
ciple forbidden in ccntro-symmctric media |23l \24\ . It is 
therefore particularly suitable to follow the re-orientation 
dynamics between an initial polar state and a final non- 
polar state. The principle of the experiment is the fol- 
lowing: the films to be studied are first brought in the 
liquid phase and a strong DC poling field E p , parallel 
to the film plane, is then switched on. This field tends 
to align the dipoles carried by the side-groups along the 
direction of the field (Fig.[3J). In this polar state, the sam- 
ple generates a high second-harmonic (SH) signal. At a 
time t = 0, we switch the poling field off and follow the 
reorientation of the dipoles toward an isotropic in-plane 
distribution of dipoles orientations through the decay of 
the second-harmonic signal in time. Other glassy sys- 
tems have already been studied with this experimental 
method [U IS S IH , most of them being poly- 
meric materials. In 25] for instance, the measurements 
probe the reorientation dynamics of a chromophore that 
is either dissolved in the hosting material or attached to 
the polymer main chain by a flexible spacer. 

In the following, we start by describing our samples 
and experimental technique (Section [HJ. We present 
then our results on the re-orientation dynamics in the 
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FIG. 3: Schematic representation of the orientation order of 
the polar branches in the glass former with applied poling field 
(left) and after removal of the field and complete relaxation 
(right). The poling field Ey is perpendicular to the opti- 
cal plane (XZ) defined by the incident and second-harmonic 
beams. 



glassy films and how it is affected by changes in the size 
of the system (Section and finally discuss these re- 
sults (Section EJ. 



II. EXPERIMENTAL 



the middle of the gap the electric field is minimum and 
equal to E p = -K, which corresponds to 1.2 -10 6 V/m. 
Around this minimum the field is constant within 1% in 
a region of width d/5 — 200 /xm in our case. This width 
is sufficiently large for our measurements, which are all 
performed with a projected laser beam area of ^100 /im 
in diameter. We assume therefore that the electric field 
is uniform in the probed surface area. The degree of po- 
lar ordering induced by the poling field is proportional 
to the applied field. Therefore a higher field is expected 
to increase the ordering and the initial SH signal before 
the field is switched off. Nevertheless we have to limit 
the applied voltage to 1.9 kV because higher fields cause 
arching between the two electrodes, which destroys the 
samples. The high-voltage power supply limited the cur- 
rent during the transient regime just after switching the 
electric field on or off to 10 mA. 

To eliminate the solvent from the deposited films, sam- 
ples were first annealed at 45 °C for 5 hours. They were 
then heated above the glass transition to 120 °C where 
they were kept for at least 15 minutes. After switch- 
ing on the external field the temperature was brought 
to the desired value and then kept constant during the 
experiment. All temperature ramps were conducted at 
0.25 Kmin -1 . To facilitate solvent removal and provide 
a clean atmosphere, the sample chamber was fluxed con- 
tinuously with dry nitrogen. 



A. Samples 



For the relaxation measurements we have used the ma- 
terial in FigH(T g = 106/110 °C) [Hill, that belongs 
to the polyalcohols family. This family is characterized 
by the presence of hydroxy functional groups that are 
known to create H-bonded networks. In the plot of the 
relaxation time versus T g /T, as originally done by An- 
gell 0, , alcohols j3^] and polyalcohols [34| are found 
to occupy the intermediate region between strong and 
fragile glass-formers. 

Thin films were obtained by spin coating from a cy- 
clopentanone solution on fused quartz substrates. The 
thickness (100-1200 A) and quality of these films were 
characterized by x-ray reflectivity |22t I35I ] and atomic 
force microscopy j^- It had already been observed by 
x-ray reflectivity that films of the same compound with a 
thickness between 40 and 100 A become unstable during 
the first heating and break up into droplets j3|J . 

To be able to apply an electric field parallel to the sub- 
strate, the quartz plates were first patterned with two 
planar electrodes separated by a gap d of 1 mm. The 
electrodes were obtained by first depositing a thin (10 
nm) layer of chromium, which has a high affinity for sil- 
ica and then depositing a second thick (100 nm) gold 
layer. A high DC electric field was applied between the 
two electrodes. For poling the films, we used a constant 
voltage V of 1.9 kV. The electric field produced in the 
gap is perpendicular to the electrode edges and paral- 
lel to the substrate surface. The intensity profile of this 
field depends on the distance from the edges [13 ■ In 



B. Measurements 

The reorientation of the dipoles has been followed mea- 
suring the second-harmonic signal after switching off the 
poling field E p . We have measured the second-harmonic 
signal with an in-coming light and SH signal having both 
an s-polarization (s—Y in Fig. [21 s-in-s-out) perpendic- 
ular to the incident plane (XZ in Fig. By releasing 
of the field at time t — 0, the field- induced orientational 
order disappears in the course of time by the randomiza- 
tion of the orientations of the dipoles. This leads to a 
decrease of the SH signal to zero. So the reorientation 
dynamics of the system can be followed by measuring the 
SH signal I ss . The square root of the I ss is directly pro- 
portional to the degree of polar ordering in the system 
through the effective susceptibility Xss of the film [38| . 
The relaxation functions are determined over a range of 
11 orders of magnitudes in time from 1 /is upward. The 
ability to characterize reorientation dynamics over a wide 
range of time scales is crucial to determine accurately the 
fitting parameters of the relaxation function. 

With the above measurement method, the state from 
which the relaxation starts is determined by the applied 
poling field. For the range of poling fields we have used, 
X ss (0) is proportional to the field strength [2(],|3!l,E3. ^ n 
principle, the measured relaxation times can be expected 
to depend on the degree of order of the starting state of 
the relaxation, thus on E p |4l|. We have not studied this 
dependence and have used only one value of E p . 
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For the second-harmonic generation measurements we 
have used a frequency-doubled Q-switched mode-locked 
Nd-YAG laser (Antares 76, Coherent) as light source with 
an output wavelength at 532 nm, using a s-in-s-out po- 
larizations configuration (for more details see |35j). The 
repetition rate of the Q-switched pulses was tuned be- 
tween 50 and 200 Hz depending on the time scale at 
which the relaxation process was explored. A special 
setup was developed to trigger and delay the switching 
of the external poling electric field with respect to the 
optical Q-switched pulses in order to reach time scales 
for the relaxation process down to 1 (is [35l E^ . The 
Q-switched pulse energy was at most 1 mJ, which is not 
high enough to induce molecular reorientation. We have 
also checked that the laser does not influence the dynam- 
ics of the system by heating it. All measurements were 
averaged over a few thousands repeated collections to in- 
crease the signal to noise ratio using a home-developed 
collection protocol |35j . 

C. Data analysis 

The measured values of the effective susceptibility 
Xss of the film are normalized to the value before the 
switching off of the poling field (xss(0)) and fitted using 
the phenomenological stretched exponential Kohlrausch- 
Williams- Watts (KWW) response function 0: 

Xso® = = e-W"?* (1) 

Xss{0) 

where tw is a characteristic decay time and < (3w < 1 
is the stretch exponent. This non-exponential relaxation 
function is generally associated to a distribution of pur e 
exponential relaxations spread through the system |43| . 
The system is then considered as a set of relaxators i, 
which can be visualized as nanoscopic regions within the 
glass- former, each one with a defined time constant Tj. 
The observed relaxation results from their superposition: 



TABLE I: Values of the characteristic time tw and stretch 
exponent f3w obtained by fitting the experimental relaxation 
with Eq. Q, and average relaxation (r) time calculated with 
Eq. for the different temperatures and thicknesses con- 
sidered. 

h(0) T (°C) t (sec) fi (r) (sec) 

116 nm 60 1.46 -10 6 0.26 29 -10 6 



70 nm 
41 nm 
25 nm 
17 nm 



10 nm 



by Emri and TschoeglJ0L_and often applied to inter- 
pret SHG decays [20 L l46 t |47 | . The maximum of G(r) 
corresponds approximately to Tw and its breadth is re- 
lated to < Pw < 1< The smaller (3\y the wider the 
distribution, the special case fiw = 1 corresponding to a 
single-exponential relaxation and a delta function G(r). 



III. RESULTS 
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XNo(i) = X)G(r i )e-*^ (2) 

i 

where the weight of relaxator i is given by G(ji), more 
often written as a continuous distribution function G(t) 
where r is a variable for the relaxation times and is not 
to be confused with the KWW parameter Tw- 
in the general case when the form of XNo(i) is arbi- 
trary, the two functions Xno(^) and G(r) are uniquely 
derived from one another through an operation involving 
a Laplace transform In the case when G(r) follows 
the KWW law, G(r) is directly determined by the two 
parameters tw and 0w- G(t) can then be written in 
the form of a definite integral or a series, but there is no 
analytical solution except for the cases flw — 1/2 and 
Av = l. 

The distribution function G(r) is numerically com- 
puted using a recursive computer algorithm developed 



We have examined the dependence of the reorienta- 
tion dynamics as a function of two parameters: tem- 
perature and film thickness. In the following, we first 
present the results concerning the temperature depen- 
dence of the dynamics in a thick film (as-spun thickness 
/io=116 nm). From the thermal expansivity measure- 
ments we have presented in [22J, we expect a bulk-like 
behavior for the molecular dynamics of such a thick film. 
We can therefore compare our dynamic measurements 
with the expected behavior around the glass transition 
temperature of a bulk glassy system. Next we examine 
how the dynamics at different temperatures is affected as 
the film thickness decreases. From the broadening of the 
glass transition we have reported in |22j . which is a sign 
of an increased inhomogeneity of the film, we expect a 
broadening of the distribution of the relaxation times as 
the film thickness decreases. 
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FIG. 4: Normalised second-order susceptibility vs. time for a 
116 nm-thickfilm (as-spun thickness) at (o) 60 °C (T s -46) and 
(A) 100 °C (Tg-6). Solid lines represent the best least-square 
fit using Eq. Q . 
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FIG. 5: Distribution of relaxation times calculated for a 116 
nm-thick film from the parameters given in Table Q] 60 °C 
(o), 80 °C (■), 90 °C (o), 100 °C (A) and 120 °C (•). The 
solid lines are only a guide to the eye. 



Temperature-dependence of relaxation 
dynamics in thick films 



In Fig. 2] we plot the time dependence of the effective 
non-linear susceptibility \ ss normalized to its value be- 
fore the poling field is switched off at time t = 0, for a 
116 nm-thick film. We show for clarity the results only 
for the temperatures just below (100 °C) the nominal 
bulk glass transition temperature determined by differ- 
ential scanning calorimetry (DSC) (T ff =106 °C) and deep 
into the glass phase (60 °C). From the observed decays 
it is clear that any contribution to Xss from electric-field- 
induced second-harmonic (EFISH) generation [4H 1421 l48l ] 
can be neglected. Indeed, for all temperatures, the SH 
signal does not decrease for several decades in time. A 
possible contribution of EFISH would result in an im- 
mediate jump of the signal to a lower value, as fast as 
the poling field E p drops to zero (~ 100 ns). We con- 
clude then, that the observed decays are associated solely 
to reorientation dynamics of the non-linearly polarizable 
units of the molecules. 

Data are fitted with the stretched-exponential KWW 
law given by Eq. l[T]l. The characteristic time tw and 
the stretch exponent Pw obtained at the different tem- 
peratures are given in Table [I] The values of (3w exhibit 
little temperature dependence. As they are significantly 
different from 1, it is obvious that the relaxation process 
is highly non-exponential. 

Another way to characterize the dynamic behavior of 
the system is to calculate an average reorientation time 
(t) as is generally done for polymer systems from dy- 
namic optical or dielectric measurements [l(J, l2a bM 



(r) 



exp 



t y w 

TW ) 



t w T (1//3 W ) 

dt = o ( 3 ) 

Pw 



where Y is the gamma function. For (3w = 1 (single 
exponential), (t) is equal to tw- The results for the 116 
nm-thick film are given in Table [fl For all temperatures 
(t) is rather large. Even for the highest temperature 
(T = 120°C= T g + 14°C) the average relaxation time lies 
above the value generally associated to the macroscopic 
glass transition (t)t s = 100 s Q. 

From the KWW fitting parameters tw and (3w we 
have computed the distribution functions of relaxation 
times G{t) as explained in Section III CI (Fig. [SJ ■ The 
measured relaxation decays at low temperatures do not 
decay to at the largest experimental time (Fig.^J). Cal- 
culating distribution of relaxation times using then high 
characteristic time tw evaluated with the KWW fitting 
law (Eq. J2J) introduces weights G(ri) for times r, longer 
than the measured ones. Thus the contributions to the 
distribution function at times log(r) > 6 have to be taken 
as purely a result of the computation but not directly 
related to the experimental measurements. The calcula- 
tions of G(t) confirm the trend indicated by Table [I] for 
the average relaxation time and the width of the distribu- 
tion. For the considered film thickness, the width of the 
distribution does not show a particular trend at tempera- 
tures larger than 80 °C. Only at the lowest temperature 
measured (60 °C) a slight widening of the distribution 
is visible. When the temperature increases, the peaks 
of the distribution move from long to short times by al- 
most four orders of magnitude in time in the explored 
temperature window. Note that the characteristic time 
Tw corresponds approximately to the maximum of the 
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FIG. 6: Xno vs - time for a 17 nm-thick film (as-spun) at (o) 
60 °C (T s -46) and (•) 120 °C (T g +U). Solid lines represent 
the best least-square fit using Eq. Q. 
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FIG. 7: Distribution of relaxation times for a 17 nm-thick film 
from the parameters in TableQ 60 °C (o), 90 °C (o), 100 °C 
(A), 110 °C (□), 120 °C (•) and 130 °C (V). The solid lines 
are a guide to the eye. 



distribution, while the average calculated with Eq. © 
is equal to the average calculated from the distribution 
G(r). This confirms the consistency of our data analysis. 

The wide distributions of relaxation times shown in 
Fig. [5] confirms the fact that to make an accurate eval- 
uation of the distribution G(r), experimental relaxation 
decays covering many orders of magnitude in time are 
needed. Indeed, the relaxation time distributions cover 
14 orders of magnitude in time (10 -4 — 10 10 s) for the 
studied range of temperatures. 



B. Thickness-dependence of the relaxation 
dynamics 

To investigate the effect of confinement we carried the 
same relaxation experiment as presented in Sec. lIII Al for 
films with a lower thickness. In Fig. we show the signal 
decay measured at 120 °C and 60 °C for a 17 nm-thick 
film (as-spun thickness). 

The decay measured at 120 °C presents an interesting 
feature: the SH signal exhibits a first decay in the short 
time domain. The use of the KWW equation to fit the 
data appears to be only a first approximation. The data 
are not properly described in the ms time window and a 
general double stretched exponential law would be more 
suitable. It seems that, under our experimental condi- 
tions, it is possible to discriminate different relaxation 
processes. Measuring relaxation dynamics around and 
below the bulk glass transition means in principle probing 
the a-relaxation process. Nevertheless, at high temper- 
atures, the a-process does not necessary screen all other 
faster processes. These modes observed at high temper- 
ature can have different origins. One possible interpre- 
tation is to attribute this mode to the fast /3-process, 
related to a local mobility. 

However, given the small contribution of possible fast 
processes (observed only at high T), we fitted in first 
approximation all decays with a single stretched expo- 
nential function using Eq. Q. Fitting with a double 
stretched exponential the data at high T hardly changes 
the parameters for the a-process. The resulting fitting 
parameters tw and (3w together with the calculated av- 
erage relaxation time (r) are reported in Table [I] The 
corresponding distribution of relaxation times G(r) are 
plotted in Fig. |7| 

The distribution obtained in the glass phase are 
slightly wider than the ones obtained with the 116 nm- 
thick film. They also show a clear increase of the expo- 
nential character (increase of (3w) for the highest temper- 
ature measured, which results in a sharper distribution 
of relaxation times for T=130 °C in Fig.0 Note that for 
this temperature the distribution function moves towards 
shorter times and (r) < 100 s. 

To increase the effect of confinement, we repeated the 
above experiments with an even thinner film. The only 
difference with the thicker films was that the annealing 
temperature was 110 °C instead of 120 °C, due to the 
instability of the film at higher temperatures. 

Fig. [H] shows the relaxation for the 10 nm-thick film 
at 90 °C and 110 °C. For this film thickness experimen- 
tal data are in general of lower quality with a higher 
scattering of the points. This is mainly due to the low 
value of the SH signal intensity in the oriented state at 
t < 0. Analysing ultra-thin films reduces the number of 
molecules probed by the input laser beam, reducing the 
sensitivity of the technique. The fitting parameters of the 
SH decays for four different temperatures are presented in 
Table [I] The calculated distribution of relaxation times 
are shown in Fig. |5| 
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IV. DISCUSSION AND CONCLUSION 
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FIG. 8: Xno vs - time for a 10 nm-thick film (as-spun thick- 
ness) at 90 °C (T 9 -16) (o) and °C (T s +4) (□) 110. Solid lines 
represent the best least-square fit using Eq. Q. 
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FIG. 9: Distribution of relaxation times for a 10 nm-thick film 
from the parameters in Tabled 60 °C (o), 90 °C (o), 100 °C 
(A) and 110 °C (□). The solid lines are a guide to the eye. 



As temperature increases, the peak in the distribution 
shifts uniformly towards shorter times till 100 °C. As 
for the 17 nm-thick film, f3w first decreases and then 
increases slightly as the temperature decreases. In both 
cases, the minimum in f3w is at around 100 °C (T g -6 °C). 

At this temperature we have also collected the relax- 
ation decay for 25 nm-thick, 41 nm-thick and 70 nm-thick 
films. The distribution of relaxation times G(t) of all the 
measured films at 100 °C are plotted in Fig. 1101 



We summarise our results on the dynamics of thin films 
by plotting the temperature dependence of the average 
relaxation time (r) (Fig. Illfl and the stretch exponent 
fiw (Fig- HH for all the films we have measured. 

The first important issue arising from these measure- 
ments is whether the measured reorientation dynamics of 
the non-linearly polarisablc groups is coupled with any 
glassy dynamics of the molecules as a whole. An ele- 
ment of answer can be found in the values of (3w for the 
thickest film exhibiting a bulk-like behaviour. The values 
are rather scattered but remain in the interval 0.25-0.35. 
If we take for the value at T g the average of the mea- 
sured values, we obtain (3 w (T g ) ~ 0.30, which is close to 
the value for bulk poly (propylene glycol). This can be 
expected from the fact that our compound carries four 
hydroxy groups per molecules and should belong to the 
most fragile hydrogen-bonded glass-formers. Another el- 
ement of answer can be found in the Arrhenius plot of (r) 
as a function of the inverse temperature (Fig. Illj l where 
(t) increases exponentially with 1/T as expected for a 
glassy material. 

We thus conclude that the reorientation dynamics we 
observe, corresponds to the one expected for a glass- 
former. The reorientation of the non-linearly polaris- 
able groups is probably strongly coupled to that of the 
rest of the molecules, since the flexible spacer between 
the central rigid group and the polarisable side-groups 
contains only three atoms (see Fig |2J) . The fact that 
the reorientation dynamics of side-groups is coupled to 
the a-relaxation process has also been demonstrated in 
a polymeric system with a flexible spacer of four atoms 
[4fij . The same coupling has also been demonstrated for 
chromophores dissolved in a polymer without any cova- 
lent bonding with the polymer chains |2q . However this 
could be due to strong steric interactions between the 
chromophores and the chains, which are less likely to 
take place in our low-molecular weight glass- former. 

We also observe however deviations of the dynamics 
behavior from the expectations on two points. One is 
that even at temperatures above the bulk glass transi- 
tion temperature T g (as measured by DSC), the average 
relaxation time (r) is longer than the one corresponding 
to the standard definition of T g : {t)t 9 — 100 s. In we 
saw that there is no evident variation of T g under confine- 
ment. So the high values of (r) cannot be explained by 
an increase of T g . The origin of the long relaxation times 
might be related to the fact that the non-polarisable side- 
groups are polar mesogenic units. Ordinary non-glassy 
liquid crystals based on the same cyanobiphenyl units ex- 
hibit also a very slow reorientation dynamics in ultra-thin 
films (thickness below 10 nm) [43, 03> 03 • The reason 
for this slow dynamics is not clear. One possible ex- 
planation could be the electrostatic correlation between 
the dipoles carried by the cyanobiphenyl groups |40| , al- 
though no direct evidence for such a correlated dynamics 
has been found yet. 
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FIG. 10: Distribution of relaxation times at 100 °C (T s -6) for 
a 10 nm (□), 17 nm (+), 25 nm (o), 41 nm (x), 70 nm (o) 
and a 116 nm (A) thick films. The solid lines are a guide to 
the eye. 
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FIG. 11: Average relaxation time (r) vs. inverse temperature 
for 10 nm (□), 17 nm (+), 25nm (o), 41 nm (x), 70 nm (o) 
and 116 nm (A) thick films. 



The second unexpected feature in the observed dynam- 
ics (in particular in thick films) is the fact that the time 
dependence of the observed relaxations remains strongly 
stretched at temperatures above the bulk glass transi- 
tion temperature. We could expect an increase of the 
stretch exponent 0\y indicating a narrowing of the dis- 
tribution as temperature increases. The fact that we do 
not observe this narrowing might simply be due to the 
limited range of temperature that we can explore above 
T g . When we have managed to performed experiments at 
slightly high temperatures (130 °C with the 17 nm-thick 
film), we indeed observed an increase of w ■ Similar re- 
sults were found on polymer systems [26l wK, l54| , where 
0w does not vary significantly in the glassy phase, while 
for temperature above T g , 0w was increasing with tem- 
perature. 

Our experimental results allow a discussion on the 
main purpose of our study, namely the effect of confine- 
ment on the dynamics of glass-formers. Let us start with 
the average relaxation time (r) (Fig. 1111) . At high and 
low temperature, (r) hardly depends on the film thick- 
ness. At intermediate temperatures around T g , there is a 
large scatter without a specific variation with film thick- 
ness. This is likely due to the limited relevance of (r) 
when the distribution of relaxation times G(r) becomes 
very wide. Variations in the actual shape of the distribu- 
tion can have a strong influence on (r) while the location 
of the maximum of G(t) hardly moves. This can be seen 
in the plot of the distribution of relaxation times for dif- 
ferent film thickness at T = 100 °C (Fig. EH). Except for 
the thickness of 17 nm, maxima are hardly shifted with 
respect to one another (less than one order of magnitude), 
but the asymmetric widening of G(r) for the thin films 
leads to differences in (r) of two orders of magnitude. So 



the average relaxation time (r) is not the most relevant 
parameter to describe the variation of the distribution of 
relaxation times with film thickness. One should rather 
look at the width of the distribution, which is related to 
the stretch exponent 0w ■ The summary of the results in 
Fig. 1121 show that we can distinguish three temperature 
regions exhibiting different behaviours: 

• high temperatures (T > 120 C): although the 
number of measurements we could perform in this 
range is small, it is clear that thickness has little 
influence on the dynamics. This can be understood 
from the fact that at 120 °C, in the liquid phase, the 
dynamics can only be little affected by size effects. 
Eventually the dynamics will be modified at small 
thickness (below 10 nm) by the presence of the free 
surface and film-substrate interface. The effect is 
quite likely not large enough to be observed with 
our experiments. 

• intermediate temperatures (70 °C < T < 120 ° 

C): in this temperature range, there is an influence 
of the film thickness on the dynamics. First of all 
the stretch exponent 0w of the thinnest films (10 
nm and 17 nm) significantly decreases as temper- 
ature decreases and reaches a minimum at around 
100 °C and then slightly increases again (Fig. 112(1. 
The decrease of 0w indicates a widening of the dis- 
tribution of the relaxation times. This widening is 
also observed for intermediate thicknesses at T = 
100 °C, but the effect is smaller and decreases as 
the film thickness increases. Since this widening is 
not observed in the thickest film (116 nm-thick) we 
can conclude that it is due to an increase of the in- 
homogeneity of the films under confinement. This 
is in agreement with our observation of a broad- 
ening of the glass transition as the film thickness 
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FIG. 12: Stretch exponent (3w vs. temperature for 10 nm (□), 
17 nm (+), 25nm (o), 41 nm (x), 70 nm (o) and 116 nm (A) 
thick films. 



decreases. Note also that the temperature range in 
which f3\y is decreased in the thinnest films, with 
respect to the value in thick films, corresponds to 
the temperature range in which the glass transi- 
tion takes place in the thinnest films studied with 
x-ray reflectivity . This is indeed the tempera- 
ture range in which one expects the inhomogeneity 
to be the largest, with part of the film in the liq- 
uid state and the rest in the glassy state. Since 
the average relaxation time varies very rapidly in 
the vicinity of the glass transition, the difference 
between the dynamics in the different parts of the 
system is then the largest. 

• low temperatures (T < 60 °C): the difference in 
dynamics between films at different thicknesses di- 
minishes as temperature decreases below 100 °C 
and becomes unmeasurable at approximately 60 
°C. Although we have not performed measurements 
at lower temperatures because the relaxation of the 
system becomes extremely long, we can expect that 
this independence on temperature also occurs be- 
low 60 °C. 

Following the argument in the discussion of the in- 
termediate temperatures, we can say that 60 °C 
is the temperature at which the glass transition is 
completed in the whole film (for the thinnest films 
we have measured). The inhomogeneity of the dy- 
namics induced by confinement becomes then equal 
to or smaller than the intrinsic inhomogeneity of 
the dynamics in bulk and thick films. 

Our finding that flw becomes independent of film 
thickness at low temperatures is is contrast with the 
results of Hall et al. [2(j on polymer films. They 
found little influence of temperature on the differ- 
ence in value of fiw for thin and thick films. This 



might however be due to the fact that they only 
performed measurements for temperatures above 
T g — 10 °C, while we measured down to T g — 46 
°C. 

In the above discussion, we have not considered the 
role of the cooperative length £ and its temperature de- 
pendence. From the model presented in |55| . we know 
that the existence of £ introduces an inhomogeneity of 
the system for film thicknesses close to £. The origi- 
nal theory of Adam and Gibbs of cooperative motion [6| 
predicts a divergence of £ at a temperature below T g . 
This implies that £ should increase as temperature de- 
creases. This means that as temperature decreases, the 
effect of confinement on the heterogeneity of the dynam- 
ics should start appearing at larger thicknesses and, for 
a given thickness, the inhomogeneity should increase as 
temperature decreases. This is however not what we have 
observed. One explanation for this is that £ does not in- 
creases as temperature decreases, or at least not much. 
This has actually been observed in colloidal suspension 
using confocal microscopy 0, |5|J- In such a system, 
the glass transition is obtained by increasing the den- 
sity of particles. The size of fast relaxing clusters (which 
can be compared to cooperatively rearranging regions) 
was found to remain constant when the density increases 
beyond the glass transition. Unfortunately, no data is 
available on the temperature dependence of £ below T g 
in molecular systems. 

In the simple model we presented in [liij . making the 
cooperative length £ independent of temperature, makes 
also the relaxation time in bulk and the degree of inho- 
mogeneity in the film independent of temperature. We 
could introduce a variation of relaxation time in bulk 
with temperature in a phenomenological way by impos- 
ing the appropriate temperature dependence of the en- 
ergy barrier per molecule that needs to be overcome for 
allowing for a relaxation process ^3[- There is however 
no obvious way of obtaining a dependence of the distri- 
bution of relaxation times in the films as a function of 
temperature in agreement with our experimental results. 
So the simple model based on the idea that motion is co- 
operative over a certain length scale is not able to explain 
all the dynamic feature of thin glassy films, although it 
explains qualitatively some of them. 

In conclusion, we have investigated how the distribu- 
tion of relaxation times in a low-molecular- weight glass- 
former is affected by confining the sample in a thin film. 
Our main finding is that the effect of confinement on the 
dynamics of the system is the largest around the glass 
transition temperature in bulk and that this effect van- 
ishes as temperature is further lowered deep in the glass 
phase. 
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